Electron Temperature of Ultracold Plasmas 
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We study the evolution of ultracold plasmas by measuring the electron temperature. Shortly 
after plasma formation, competition between heating and cooling mechanisms drives the electron 
temperature to a value within a narrow range regardless of the initial energy imparted to the 
electrons. In agreement with theory predictions, plasmas exhibit values of the Coulomb coupling 
parameter F less than 1. 
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Knowledge of the electron temperature is critical to 
understanding the behavior of ultracold plasmas formed 
by the photoionization of laser-cooled atoms Q- Un- 
til now, this important parameter had not been directly 
measured. Simple consideration of the photoionization 
process (i.e. considering only single atoms in isolation) 
implies that the initial kinetic energy of the electrons in 
the plasma is proportional to the energy of the photoion- 
ization photon in excess of the ionization limit. In this 
case, only the linewidth of the photoionization laser limits 
the lowest energy imparted to the electrons, and in prin- 
ciple plasmas with electron temperatures below IK could 
be created. The relationship between the actual electron 
temperature and temperatures obtained from this simple 
consideration remains an open question. 

Several strong heating and cooling processes in the 
plasma can have a radical effect on the temperature 
of the system. Continuum lowering 0, correlation- 
induced heating 0, 0| , three-body recombination to Ry- 
dberg atoms and the evolution (deexcitation) of the Ry- 
dberg atoms |E S S S ^| are all predicted to result 
in heating of the electrons that increases with increasing 
density. Conversely, the electrons possess thermal pres- 
sure that causes the plasma to expand and that in 
turn induces strong adiabatic cooling. Evaporative cool- 
ing is certainly present as well. All of these processes can 
act within the first few microseconds after the plasma is 
created, and so the evolution of the electron temperature 
will be determined by the balance of these (and perhaps 
other) heating and cooling processes. 

The determination of the electron temperature will be 
useful in interpreting the results of future experiments 
with ultracold plasmas. Measurements of the electron 
temperature test molecular dynamics simulations' 0, Q 
predictions of this fundamental plasma parameter. Be- 
cause electron-atom collisions such as Rydberg recom- 
bination and superelastic scattering influence tempera- 
ture, atomic collision theory in plasmas can also be tested 
Pl la y. Ilil • Finally, a measurement of the electron tem- 
perature is necessary in order to determine to what ex- 
tent, if any, the electrons of these ultracold plasmas are 
in the strongly coupled (i.e. correlated) regime [l^ . 

The ratio of Coulomb energy to thermal energy, F, 



expresses the importance of correlations in the plasma. 
This parameter is defined (for the electrons) as F = 

— !— where a — i^r—Y^^ is the Wigner-Seitz ra- 
dius, e is the electron charge, n is the average plasma 
density, and T is the electron temperature. For F > 1 
(the strongly coupled regime), spatial correlations in the 
plasma develop and phase transitions are possible |l2l |. 
It is predicted 0, S S 8, 9] that various heating ef- 
fects in the plasma will each limit F to be less than 1, 
so measuring F also tests the validity of these theoretical 
treatments. In general, strongly coupled neutral plas- 
mas (distinct from strongly coupled non-neutral plasmas 
[l^l Im ll^L llfij |) are rare and difhcult to produce in the 
laboratory; ultracold plasmas could potentially be cold 
enough to reach the strongly coupled regime. 

In this work we measure the electron temperature a 
short time after the plasma formation. The creation of 
ultracold plasmas in our apparatus has been previously 
described in Ref. 0. Briefly, we use a Magneto-Optic 
Trap (MOT) to collect and cool 4x10^ metastable xenon 
atoms to a temperature of about 20/iK. The spatial den- 
sity distribution is roughly gaussian with a typical rms 
radius a ~ 250 /im. We produce the plasma in a 10 ns 
two-photon excitation of up to 30% of the initial sample, 
with the number TV of photoionized atoms and the initial 
electron energy A_E controlled by the intensity and fre- 
quency of the photoionizing laser, respectively 0. We 
apply small electric fields (~ 0.3 — 4 V/m) to the plasma 
during its evolution by controlling voltages on wire mesh 
grids (>95% transmission) positioned above and below 
the plasma region. Electrons that leave the plasma are 
directed to a microchannel plate (MCP) detector. We 
determine the density and size of the plasma by measur- 
ing its radio-frequency (RF) response as in [ill, with a 
slight modification to the calibration used there. In light 
of recent theoretical work 0, , we now base our size 
calibration on the measured asymptotic expansion veloc- 
ity of high A£: {AE/kB > AOOK), low (N < 5x10^) 
plasmas |l9j . 

After the photoionization pulse, some electrons rapidly 
escape as a result of the kinetic energy imparted dur- 
ing photoionization. This results in an excess of ions as 
compared to electrons in the plasma region and a space 
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FIG. 1: Self-consistently calculated electron potentials show 
as a function of position along the applied electric field dire^ 
tion, for ^=1.5 V/m (solid line) and 3 V/m (dashed lim 
applied electric fields. The relevant plasma parameters ai 
N,on = 10^ Neiectron = 0.95N,on, <T = 0.53mm, aud T=40i 
The value of the parameter Ecut{F) used in the calculation ( 
the fraction of electrons spilled is determined from the sadd 
point as indicated. 



charge develops. This space charge in turn confines tl 
remaining electrons in a potential well Elastic coll 
sions rapidly redistribute the remaining electrons' energ 
into a thermal distribution |20l |. Since the ions in tl 
plasma are unconfincd, the plasma expands outward i 
response to the thermal pressure of the electrons [llj . 

Due to the low temperatures and small size of tl 
plasma, the standard methods used to measure plasrr 
temperatures are not applicable [iJl- We developed 
method in which we probe the energy distribution of tl 
electrons to measure T. After a desired amount of ev( 
lution time, an electric field is turned on. This field sut 
denly lowers the lip of the potential well confining tl 
electrons, and energetic electrons become unconfincd (sc 
Fig. These electrons then spill from the plasma an 
strike the MCP. Figs. |21 and |21 show a typical spill signi 
and a plot of the fraction of electrons spilled as a fun( 
tion of the electric field amplitude F, respectively. I 
order to avoid having to dynamically model the chang- 
ing confinement during the spill, the fraction spilled / is 
kept small enough so that the ratio N electron /N ion is not 
altered substantially. 

A model is necessary in order to extract the temper- 
ature from the measured spilling curves. We developed 
a method to solve for the electron density distribution 
given the ion density, electron number, and temperature. 
Since the electrons' confining potential is influenced by 
the electrons themselves, their density distribution is cal- 
culated self-consistently, starting from an approximate 
distribution and iteratively improving that distribution 
until it converges to a self-consistent result. While this 
self-consistent calculation allowed us to calculate / as a 
function of F, we found that / calculated this way was 
highly sensitive to the details of the ion density distri- 
bution. We did not want the measured temperature to 
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FIG. 2: MCP signal (solid line) showing a typical plasma 
evolution and spilling signal (region indicated in oval). The 
inset shows the spilling part of the curve in more detail. The 
magnitude of the applied electric field is also plotted (dashed 
line). 

be contingent on a simulation of the motion of the ions 
in the plasma (especially since the density of atoms in 
our MOT is not perfectly spherically symmetric or gaus- 
sian). Therefore, we switched to a more robust method 
that does not depend explicitly on the ion density distri- 
bution. This method relies on the fact that at positions 
far enough away from the center of the plasma, the elec- 
tron potential is U{r) — _ (^'°''°~-^^^^'^*'-°"°)'^ 1^ regard- 
less of the ion and electron density distributions. The 
fraction spilled as a function of F is then modeled using 
the following expression: 

/ = /3 + a/ Z?(6)exp(--— )de (1) 

where D{e) is the density of states calculated for the 
plasma's electron potential in its asymptotic limit, A and 
[3 are fit parameters corresponding to the chemical poten- 
tial and a background subtraction, Ecut {F) is determined 
by the magnitude of F as shown in Fig. ^ and Ecut{Fo) 
is determined in the same way for i^o, which is the bias 
electric field present before the spilling field is turned on. 
As long as the saddle point associated with Ecut is far 
enough from the center of the plasma, Eq. (1) will be a 
good approximation of the number of electrons spilled in 
response to lowering the potential barrier. 

We tested the validity of this model in several ways. 
First, we compared values of T derived from fits to 
Eq. (1) with T determined from the more sophisticated 
self-consistent model for a variety of different ion den- 
sity distributions. The values of T from the two different 
methods matched to better than 20% for the plasma pa- 
rameters examined in this work. The form of D(e) was 
varied and found to have little effect on the fit tempera- 
ture [2^ . Finally, experiments were performed where RF 
pulses were applied to impart a variable amount of heat. 
The temperatures determined using Eq. (1) shortly after 
these pulses scaled linearly with applied RF power 23] . 
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FIG. 3: Fraction of electrons spilled (/) vs. Ecut for a typie 
set of data (filled circles). Fits using the model based on E 
(1) are shown for the temperatures that are indicated, whe 
the 35K curve is the best fit curve and the other two a 
shown for comparison. Note that the value of T is determin 
primarily from the curvature of / vs. Ecut- 



In order to accurately calculate / in both the se 
consistent model and in the model represented by E 
(1), we must include the screening of the external field 1 
the plasma. This screening is approximated by the use 
a "screening radius," within which the external electi 
field does not penetrate. The resulting induced dipc 
field is then included in the calculation of Ecut- We e 
timate values for the screening radius by calculating t' 
first-order correction, due to the presence of the electi 
field, to the electron density distribution self-consistent 
determined in the absence of the field. 

Temperature measurements were taken over a ran 
of ion numbers (0.1-2x10^), at several different values 
AE, and at different plasma evolution times. The plasn 
evolution times that can be studied are constrained to 
from ^3-8 fis for higher AE and from '^4-15 /is for low 
AE- At early times, the electrons in the prompt pe; 
distort the spilling signal. At later times, the plasn 
expands to sizes large enough so that the saddle poir 
induced by the spilling field are inside the ^ part of t' 
potential, violating a central assumption of our moc 
and thus rendering it inapplicable. 

Figure 01 shows the measured values of T (Fig. 
and r (Fig. E^b)) vs. plasma evolution time for vario 
values of AE. Two features in Fig. ^are immediately a 
parent. First, cooling is observed as the plasma expanc 
as is expected from the dynamics of the expansion ^] 
Second, the temperatures observed for the lower values 
of AE are greater than AE itself, clearly indicating the 
importance of heating in the early stages of the plasma. 
Indeed, the range of temperatures observed falls in a re- 
markably narrow band given the factor of 30 range in 
the value of AE. As the number of ions is varied, the 
same qualitative features are present. The values of T 
decrease as N is reduced, but in such a way that F does 
not increase substantially. For the most favorable con- 
ditions (late times, AE/ks = 30K-100K, low iV) values 
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FIG. 4: (a) Measured temperatures vs. plasma evolution 
time for different values of AE/ks as indicated in the leg- 
end. These data have Nion ~ lO'', 0.4mm < cr < 0.6mm for 
AE/kB < WOK and 0.4mm < a < 0.7mm for AE/ks > 
WOK. This size range was chosen because it corresponds 
to the plasma evolution time just after the initial electrons 
have escaped from the plasma. Times 0.5-2 lis earlier (de- 
pending on A_b) and a few microseconds later were also mea- 
sured. The error bars show only the statistical uncertainty, 
(b) Corresponding F for the data shown in Fig. 4a. The 
AE/ks ~ WOK data do not continue to increase but rather 
level out at later evolution times. 



of r in the range of 0.1-0.15 were consistently achieved, 
and thus these plasmas are not in the strongly coupled 
regime at these times. It is still possible that F increases 
above 1 at much later times where we cannot apply our 
measurement methods. However, we do not see evidence 
for rapidly increasing F vs. time at the later evolution 
times we can observe. 

Several systematic uncertainties have to be taken into 
account. The dominant uncertainties are in the calibra- 
tion of the ion number (35%), in the appropriate value of 
the screening radius, and in the calibration of the plasma 
density. When combined in quadrature, these uncertain- 
ties produce a 70% overall systematic uncertainty in the 
measurement of T. 

Several assumptions are also made in the analysis: 
that the highest-energy electrons are in thermal equilib- 
rium with the rest of the electrons; that / is propor- 
tional to the number of electrons with energy greater 
than Ecut', and that the implicit truncation in the ther- 
mal distribution due to the finite potential depth does 
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not affect its Maxwell-Boltzmann character at the ener- 
gies corresponding to the values of Ecut used. By exam- 
ining spilhng curves and observing the flux of escaping 
electrons after the spilling peak, we determined that the 
plasma evolution times at which we took data were four 
times or more longer than the characteristic time that it 
took for elastic collisions to fill the spilled energy levels, 
supporting the first assumption. The last two assump- 
tions are consistent with our observations based on the 
comparison of spilling data with different ranges of Ecut 
(i. e. cutting more or less deeply for the same experimen- 
tal parameters), in which we observed no significant shift 
of T as a function of the range of Ecut- 

The observation that F is limited to values less than 
1 is consistent with theory predictions 0> S S S IS ■ 
Given the previously measured importance of Rydberg 
atom formation and evolution jQ] and the range of plasma 
evolution times studied in this work, comparing our re- 
sults to the r ~ 0.2 prediction of Refs. is more di- 
rect than comparison to the other predictions cited. For 
the relatively early plasma evolution times studied in this 
work, we measured values of T that are significantly less 
than predicted by Refs. and indeed are less than 
what would be naively implied by the other predictions 
13, S 0, 0] as well. Measurement of the full time evolu- 
tion of and T and F remains a challenge that will require 
a new measurement technique. In the future, it is pos- 
sible that measurements of the Rydberg distribution as 
a function of plasma evolution time can be used to mea- 
sure temperatures over a larger fraction of the plasma 
lifetime. 

In conclusion, we have developed a method of mea- 
suring the temperature of ultracold plasmas shortly after 
the plasmas are created. This is done by measuring the 
response of the plasma to electric fields and thereby ob- 
taining the energy distribution of the highest-energy elec- 
trons. Significant cooling was observed when the initial 
kinetic energy of the electrons was large, and significant 
heating was observed when low initial kinetic energies 
were imparted, driving the temperatures into a relatively 
narrow range. The electron components of the ultracold 
plasma were never observed to be in the strongly coupled 
regime, confirming predictions by theory. The colder ions 
may yet be in the strongly coupled regime, a possibility 
currently being studied by another group [2J|. 

The authors acknowledge useful discussions with Fran- 
cis Robicheaux during the course of this work. Robert 
Fletcher provided assistance in the collecting and analyz- 
ing the data. Scott Bergeson, Tom Killian, and Simone 
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iments. This work was funded in part by the ONR. 
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In the absence of a concrete prediction for the scaling 
of plasma temperature with applied RF amplitude, the 
linear scaling of temperature with pulse amplitude is rea- 
sonable, but not explicitly predicted. 
C. E. Simien et al, physics/0310017. 
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